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Abstract 

We study the spectral distributions of charmed meson with J p = quantum num- 
bers in cold nuclear matter applying a self-consistent and covariant many-body ap- 
proach established previously for the nuclear dynamics of kaons. At leading orders 
the computation requires as input the free-space two-body scattering amplitudes 
only. Our results are based on the s-wave meson-nucleon amplitudes obtained re- 
cently in terms of a coupled-channel approach. The amplitudes are characterized 
by the presence of many resonances in part so far not observed. This gives rise to 
an intriguing dynamics of charmed mesons in nuclear matter. At nuclear saturation 
density we predict a pronounced two-mode structure of the D + mesons with a main 
branch pushed up by about 32 MeV. The lower branch reflects the coupling to two 
resonance-hole states that are almost degenerate. For the D~ we obtain a single 
mode pushed up by about 18 MeV relative to the vacuum mode. Most spectacular 
are the results for the meson. The presence of an exotic resonance-hole state 
gives rise to a rather broad and strongly momentum dependent spectral distribution. 



1 Introduction 



It has been suggested to explore the possible mass modifications of D mesons 
with the FAIR project at GSI [1]. It is planned to study such effects by detect- 
ing D mesons in heavy-ion reactions but also in antiproton-nucleus reactions. 
The two approaches are complementary, with the second focusing on the prop- 
erties of D mesons in cold nuclear matter. Thus it is of importance to develop 
a thorough theoretical understanding of how the D mesons may change their 
properties in nuclear matter. 

At present there are only few works published on the properties of non-strange 
D mesons [2,3,4,5]. To the best of our knowledge, there is no attempt to predict 
the properties of the strange Df mesons in nuclear matter so far. Except for 
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the work by Tolos et al. [5] mean field parameterizations were considered. 
The model by Tsushima et al. [2] arrived at attractive mass shifts for the D + 
and D" mesons of about 60 MeV at saturation density. A different mean field 
ansatz by Sibirtsev at al. [3] that considers also vector mean field contributions 
suggests an attractive mass shift of about 140 MeV for the D + meson, but a 
small repulsive mass shift of about 20 MeV for the D~ meson. A QCD sum 
rule analysis of Hayashigaki [4] predicted an attractive mass shift of about 50 
MeV for the D + meson. These developments resemble to some extent the first 
attempts to predict the properties of kaons and antikaons in nuclear matter 
(see e.g. [6,7,8]), which also assumed a mean field type behavior for the mass 
shifts in nuclear matter. However, by now it is well established that for the 
antikaon such an ansatz is not valid due to complicated many-body dynamics 
induced by the presence of resonance-hole states (see e.g. [9,10]). Thus it may 
be crucial to study also the properties of D mesons in a more microscopic 
manner. 

An important step towards a more realistic approach to the D + meson in 
nuclear matter was taken by Tolos et al. [5], who did not insist on a mean 
field approach exploring the possible influence of the A c (2594) resonance on 
the nuclear dynamics of the D + meson. The key ingredient of a realistic de- 
scription are the D + N — > D + N scattering amplitudes, that determine the 
D + meson self energy at least for dilute nuclear matter [13,8]. In turn the 
results of [5] depend decisively on their particular form of the coupled-channel 
amplitudes that took into account the 7rS c , DN, r/A c channels for I = and 
7rA c , 7rS c , DN, r]T, c channels for 1 = 1. The A c (2594) resonance which carries 
J F = \ quantum numbers was dynamically generated in [5] as suggested first 
in [11]. The interaction strengths were derived relying on a SU(3) symmetry for 
the u,d and c quarks. As a consequence of their amplitudes an attractive mass 
shift of about 10 MeV together with a significant broadening of the D + is pre- 
dicted at nuclear saturation density was predicted. Recently a more complete 
computation [12] that incorporated the additional channels KE. C , KS' C , D S A 
for 1 = and KS C , KS' C , -D S S for 1=1 arrived at scattering amplitudes 
that differ significantly from those of [5]. The interaction was saturated by a 
t-channel vector meson exchange. Comparing the different interactions [5,12] 
the work by Tolos et al. severely overestimates the charm-exchange channels. 
In [12] those channels are suppressed by a kinematical factor m 2 p /m 2 D ~ 0.2. 
As an important new result of [12] the 1 = amplitude shows two resonance 
structures around 2.6 GeV. A narrow state that couples dominantly to the 
DN, D S A channels and a broader state that is interpreted as a chiral excita- 
tion of the open-charm sextet | + ground states [11]. The narrow state should 
be identified with the observed state A c (2594). The broader one, which couples 
strongly to the 7rS c channels, awaits experimental confirmation and couples 
very weakly to the DN channel. It is the analogue of the A(1405), which is a 
chiral excitation of the baryon octet | + ground states (see e.g. [14,15,16,17]). 
In contrast to [5] the isospin one amplitude of [12] reflects a narrow resonance 
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of mass 2.62 GeV which couples dominantly to the DN, _D S £ channels. This 
additional resonance will affect the properties of D + mesons in nuclear matter 
significantly. 

The purpose of this letter is threefold. First we derive the properties of the D + 
meson in nuclear matter based on the improved understanding of the D + N 
scattering amplitudes of [12]. Second, we present the mass shift of the D as 
predicted by the scattering amplitudes of [12]. Thirdly, we present for the first 
time predictions on how the strange Df mesons change their properties in a 
dense nuclear environment. Again the scattering amplitudes as obtained in [12] 
are used. The s-wave DfN — > DfN scattering amplitudes are characterized 
by exotic resonances at 2.89 GeV and 2.78 GeV. The computations are based 
on the self consistent and covariant many-body approach developed in [9,10]. 
It is important to perform such computations in a self consistent manner since 
the feedback of an altered meson spectral function on the resonance structure 
proved to be a decisive many-body affect [9,10,18]. 



2 Self consistent and covariant nuclear dynamics for charmed mesons 

In this section we recall the self consistent and relativistic many-body frame- 
work required for the evaluation of the pseudo-scalar meson propagation in 
nuclear matter [10]. The key ingredient are the vacuum on-shell meson-nucleon 
scattering amplitudes 

(Di(q) N(p) | T \D\q) N(p)} = (2tt) 4 5\q +p-q-p) 

xu(p)T% N _ ¥DN (q,p;q,p)u(p) , (1) 

where <5 4 (..) guarantees energy-momentum conservation and u(p) is the nu- 
cleon isospin-doublet spinor. Here is D = (D + , Do),(Do,DJ),Df or D~ de- 
pending on which system we study. The scattering amplitudes are decomposed 
into their isospin subsystems. The reactions DfN — > DfN carry 1 = 1/2 
inheriting the isospin of the nucleons. The amplitudes involving the isospin 
doublets (-D+, -Do) or (-Do, D ) are decomposed into isospin zero and one com- 
ponents 

n/=o) = \ 1 + (*T f ) > P SU = \ ( 3 513 1 " (*T f ) » ( 2 ) 

where q,p, q,p are the initial and final meson and nucleon 4-momenta and 
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w = p + q = p + q, k = \(p-q), k = \{p-q). 



(3) 



The in-medium meson-nucleon scattering amplitude is determined by the 
coupled-channel Bethe-Salpeter equation: 



r drl 

T(k, k; w) = K,(k, k;w) + J J^y £(^> ^ w ) ^(^ w ) ^(^ ^ w ) > ( 4 ) 

where the in-medium scattering amplitude T(k,k;w,u) and the two-particle 
propagator Q(l;w,u) depend on the 4- velocity characterizing the nuclear 
matter frame. For nuclear matter moving with a velocity v one has 



u„ = I - . = ] , u 2 = l. (5) 

We emphasize that (4) is properly defined from a Feynman diagrammatic 
point of view even in the case where the in-medium scattering process is no 
longer well defined due to a broad meson spectral function. 

To be consistent with the dynamics derived in [12] all channels must be con- 
sidered in (4) that were incorporated in [12]. In this work we exclusively study 
the effect of an in-medium modified two-particle propagator, i.e. we identify 
the in-medium interaction kernel K, in (4) with its vacuum limit to be taken 
from [12]. We restrict ourselves further to the in- medium change of the D 
mesons only, i.e. all two-particle propagators in (4) that do not involve a D 
meson are taken to be of the form used in [12]. We write 



AS(p,u) = 2-niQ(p-u) 5(p 2 - m%) + m N ) Q(k 2 F + m 2 N - (u ■ p) 2 ) , 
S(p,u) = S(p) + AS(p,u), P(g.") = g2 _ TO /_ n(g>ti) , 
G DN (l;w,u) — —iS(\w + l,u)V(\w — l,u) , (6) 
where the Fermi momentum &p parameterizes the nucleon density p with 



In the rest frame of the bulk with = (1,0) one recovers with (7) the 
standard result p = 2kp/(37r 2 ). In this work we also refrain from including 
nucleonic correlation effects. The meson self energy U(q, u) is evaluated self 
consistently in terms of the in-medium scattering amplitudes 
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— — i AS (p, u)T(±(p-q),±(p-q);p + q,u), 



(2tt)< 

r - 1 r (7=0) 4. 3 r (/=1) or r - r (/= 5 } 

J - 4 DN^DN "I" /, DN—>DN > D^N~*D^N * ^ 



In order to solve the self consistent set of equations (4,6,8) it is convenient 
to rewrite the coupled-channel system. Given our assumptions the coupled- 
channel problem reduces without any further restrictions to a single channel 
problem: 

r {i) _ T (i) T {i) ( c r \ T {i) /q\ 

1 DN^DN — 1 DN^DN ~t~ 1 DN-+DN ' [yDN ^DNJ ' 1 DN-^DN ? l y J 

where we use a compact notation. The self consistent set of equations (6,8,9) 
is now completely determined by the vacuum amplitudes T^ N _ >DN . Rather 
than specifying the coupled-channel interaction kernel of [12] it is sufficient to 
recall the s-wave scattering matrix describing the DN — > DN process. The 
one of [12] has the separable form 

r«(«,) = i + m(7) (^), 



1 ' ,(/) 



^(v/i)e 2 ^)_i) , 



2ip DA r 

Pdtv _ 1 ^ ("tjv + "flf j ^ (rag - m o) 2 j ^ 

where we recall the parametrization of the amplitudes in terms of phase shift 
5 and inelasticity rj. 

Using (10) as input in (9) the in-medium amplitude takes the form 



M {I) (w,u) = [l-M'^Vs) AJ(w,u)] 1 M<- I \yfi), (11) 

where we neglect contributions of higher partial-wave amplitudes. The p-wave 
and d-wave scattering amplitudes are not determined in [12]. The reduced 
loop function A J(w, u) acquires the generic form 



at/ x f d 4 l ,, . / I ■ w 

AJ{W > U) = J (2^ 9{l] W ' U) I mN + 77^ 
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/ x , &(k 2 F + m 2 N -(u-l) 2 ) 

g(l; w,u) = 2ne(l.u) 5(l 2 - m%) ^ 2 ' , 

v ' (w - ly - m z D — U(w - l,u) 
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l 2 — m 2 N + i e (w — l) 2 — m 2 D — U(w — I, u) 

i 1 
I 2 — m 2 N + % e (w — I) 2 — m 2 D + % e 

We observe that the loop function A J(w, u) is a scalar under Lorentz transfor- 
mation and therefore depends only on w 2 and w ■ u. Thus it can be evaluated 
in any convenient frame without loss of information. In practice we perform 
the loop integration in the rest frame of nuclear matter with — (1, 0). It is 
straightforward to perform the energy and azimuthal angle integration in (12). 
The energy integration of the last two terms in (12) is performed by closing 
the complex contour in the lower complex half plane. One picks up two con- 



tributions: first the nucleon pole leading to Iq = vm 2 N + I 2 in (12) and second 



the meson pole at typically w — l = —\]m 2 D + (I — w ) 2 . Here we neglect 
the meson pole contribution. All together one is left with a two-dimensional 
integral which must be evaluated numerically. In our simulation we restrict 
|f| < 800 MeV. 

The meson self energy follows 



d 3 p I rriN WQ — w 2 — q-w' 



n <«-"> = - 2 / wf Ur + Ep vJ ) M(w ' u) • 

M(w,u) = ^M {I=0 \w,u) + ^M {I=1 \w,u), (13) 

where we used w M = (1,0), = (ou,q) and u> M = (cu + E p ,q + p) with 
E p = [m 2 N +P 2 ) 1 / 2 . With (11), (12) and (13) a self consistent set of equa- 
tions is defined. It is solved numerically by iteration. First one determines the 
leading meson self energy IT(u;, q) by (13) with M.(w, u) = M(y / s). That leads 
via the loop functions (12) and the in-medium Bethe-Salpeter equation (11) 
to medium modified scattering amplitudes A4(w ,w). The latter are used to 
determine the meson self energy of the next iteration. This procedure typi- 
cally converges after 3 to 4 iterations. The manifest covariant form of the self 
energy and scattering amplitudes are recovered with Il(g 2 ,u;) = Il(g 2 ,g • u) 
and A4(w 2 ,wo) = A4(w 2 ,w • u) in a straightforward manner if considered as 
functions of q 2 ,u) and w 2 ,wq respectively. 
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3 Results 



We discuss the free-space scattering amplitudes derived in [12] together with 
their implications for the in-medium properties of the D mesons. The results 
are based on the self-consistent scheme [10] recalled in the previous section. 
The s-wave scattering amplitudes are shown in Fig. 1. The spectral functions 
of the D mesons are shown in Fig. 2 for meson momenta 0, 200 and 400 MeV 
and nuclear densities 0.17 fm -3 and 0.34 fm -3 . 

The amplitudes reflect the presence of various resonances. Only the D~N 
s-wave scattering process is not influenced by a resonance. In this case the 
amplitude is characterized to a large extent by the scattering length 



aglj! ~ -0.16 fm, agl# ~ -0.26 fm. (14) 

Given the isospin averaged scattering length the mass shift of a D~ meson in 
nuclear matter is fully determined by the low-density theorem [13,8]. It holds 




Fig. 1. S-wave scattering amplitudes of the D mesons off the nucleon. The amplitudes 
are taken from [12]. The normalization of the amplitudes is specified in (10). 
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Fig. 2. Spectral distributions of the D and Df mesons based on the scattering 
amplitudes of Fig. 1. Results are shown for for meson momenta 0, 200 and 400 
MeV and nuclear densities 0.17 fm -3 and 0.34 fm -3 . The self consistent many-body 
approach of [10] was applied. 

Am,- =-U— + —) (aSI + 3 aglj>) p + O (p^) 

~ 17 MeV — . (15) 
Po 

The low-density mass shift of 17 MeV is quite close to the self consistent result 
shown in Fig. 2. Self-consistency leads to a mass shift of 18 MeV at saturation 
density and 38 MeV at twice saturation density. 

The positively charged D + meson interacts with a nucleon in a more compli- 
cated manner. The s-wave scattering amplitude may be characterized roughly 
in terms of a resonance-pole term and a background term 



M(V~s)-~ 7= 7i I , - r /9 +b. (16) 

y/S - M R + I T R /2 

The isospin-zero amplitude couples strongly to the A c (2594) resonance, the 
isospin-one amplitude to a so far unobserved S c (2620) resonance: 
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a D+ N ~ -°- 43 fm > I^S?594) I - 6-6 , 

ag:^-0.41fm, 1^20)1^5.8. (17) 



In this case the low-density theorem ceases to be useful at densities much 
smaller than saturation density. A self-consistent computation is required. 
Fig. 2 demonstrates that the spectral distribution has a two-mode structure, 
which is a consequence of important resonance-hole contributions. Note that 
the masses of the A c (2594) and S c (2620) are almost degenerate. At saturation 
density the main mode is pushed up by about 32 MeV as compared to the 
free-space meson. Nevertheless, due to the resonance-hole state one may ex- 
pect that the production of D + mesons is enhanced in heavy-ion reactions as 
compared to nucleon-nucleon collisions. Our results for the D + meson differ 
from the previous study [5] significantly. This is a consequence of the quite 
different interaction used in the two computations. In particular the work [5] 
did not predict the isospin one resonance E c (2620). The latter dominates the 
resonance-hole component in the spectral distribution of the present work. 
In Fig. 3 the spectral functions of the A c (2594) and S c (2620) resonances are 
shown at saturation density as compared to their free-space distributions. 
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Fig. 3. Imaginary part of the isospin zero (l.h.p.) and isospin one (r.h.p) 
-D + -nucleon scattering amplitude at saturation density as compared to the free-space 
case (dotted lines). The amplitudes are shown for two values of the resonance 
three-momentum w = GeV and w = 0.4 GeV. 
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We observe small attractive shifts in their mass distributions and significant 
broadening, at least once the resonances move relative to the matter bulk. 

We turn to the mesons with non-zero strangeness. The study of [11] predicts 
the existence of a coupled-channel molecule in the (DjN, DA, Z)S) system. 
The latter carries exotic quantum numbers that can not be arranged by three 
quarks only. Exotic s-wave states with C=-l were discussed first by Gignoux, 
Silvestre-Brac and Richard [19] and later by Lipkin [20]. The binding of about 
190 MeV for the N s5 (2780) state predicted in [11] awaits experimental con- 
firmation. If confirmed the D~N — > D~N s-wave scattering amplitude must 
show a prominent pole structure at subthreshold energies. The amplitude in 
Fig. 1 may be characterized with 



a D-N - -°- 33 fm , \9n '$so) I - 3.1 . (18) 

The presence of the A^ sg (2780) leads to a well separated two-mode structure 
of the D~ spectral distribution. The main mode is pushed up by less than 10 
MeV at nuclear saturation density. 

Recently, it was pointed out that coupled-channel dynamics predicts attraction 
also for the (KA C , D+N, KZ C ) system [11,12]. An exotic state A 5s (2892) about 
75 MeV below the D+N threshold is predicted. From the DfN -> D+N 
scattering amplitude of Fig. 1 we extract the values 



a DtN ~ -1.21 + i O.Olfm , 1^(2892) I - 2-8 . (19) 

The in-medium spectral distribution of the derived in the self consistent 
many body approach is most striking. The two modes expected from the 
possible existence of the iV gs (2892) state are almost merged into one broad 
structure, in particular at intermediate meson momenta 400 MeV. The results 
are quite analogous to the spectral distribution of the K~ where the A (1405) 
nucleon-hole state gives rise to a broad distribution [9,10]. Clearly this result is 
an immediate consequence of the small binding energy of the A 5s (2892) state. 

We close with the quest for more detailed studies of the D-meson nucleon 
scattering processes to more reliably predict the in-medium properties of the 
D mesons. 
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